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Abstract 

We present Bc-meson two-point sum rules at next-to-leading order in NRQCD approximation. 
Analitycal results for perturbative spectral density and gluon condensate contribution with account 
for summed Coulomb corrections are given. Estimates of various c-quark masses together with 
couplings and masses of lowest lying i?c-meson resonances are performed. 

1 Introduction 



The study of the physics of i?c-mesonQ has already a long story. The collaborative efforts of many 
scientists around the world have led to remarkable predictions, describing spectroscopy [jl], 0], decays 
P, |, ^, and production mechanisms of this object. The theoretical work done was very helpful 
in the experimental search and discovery of this meson by CDF collaboration As it was many times 
earlier with other particles, the observation of Bc-meson in nature have not diminished the interest of 
physicists to this object. People wonder how much we can learn from this meson about the Standard 
Model (SM) of particle interactions. In this paper we explore a potential of Be - meson in extracting 
the parameters of SM lagragian, namely, the heavy quark masses. 

What concerns its spectroscopic properties, this meson stands among the families of charmonium cc 
r-| and bottominum bb : two heavy quarks move nonrelativistically, since the confinement scale, determining 
C^.the presence of light degrees of freedom (sea of gluons and quarks), is suppressed with respect to the 
P-i! heavy quark masses mg as well as the Coulomb-like exchanges result in transfers about asTriQ, which 
^ ] is again much less than the heavy quark mass. It is precisely the nonrelativistic nature of heavy quark 
'dynamics in the i?c-nieson, that offers us a possibility to use the NRQCD framework and gain more 
K> "inside on QCD dynamics of the constituent heavy quarks. Recently, within this framework the NLO 
*^ .and NNLO NRQCD sum rules were derived and analyzed for the case of T-family [|11[]. The result of 



5^ [this analysis was a precise determination of pole, running, IS, potential subtracted and kinetic 6-quark 
masses. Here we employ the same NRQCD sum rule method to determine the numerical values for 
different definitions of c-quark masses together with couplings and masses of lowest lying Bc-Taeson 
resonances. The present work is similar to the analysis performed for T-family with differences in 
particular analytical expressions for theoretical moments of two-point correlator. 

The paper is organized as follows. In section 2 we remind the reader the QCD sum rule framework for 
the two-point sum rule. Section 3 introduces the NRQCD framework and we comment on its connection 
with QCD one. Section 4 explains the strategy used for calculation of NRQCD two-point correlation 
function. In section 5 we calculate the perturbative theoretical expressions for moments of correlation 
function. Section 5 deals with the corrections to the correlator, coming from the gluon condensate 
operator. In section 7 we present our results on the different c-quark mass definitions, Sc-meson mass 
and coupling constant. The detail discussion on optimization methods is also given. And finally, section 
8 contains our summary. 



"'^For review see |10 
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2 Two-point sum rules 



In the QCD sum rules approach [jT2|, |T3| meson bound states are described by local interpolating currents 



of the form J = qiTq2, where F is a suitable Dirac structure to account for meson quantum numbers. 
Thus, the studied object in QCD sum rules is a coupling of meson under consideration to corresponding 
current. In the case of Be meson this coupling is defined by the following equation 



{0\bi-f5c\B,) 



nil) + 

The estimates of i?c-meson structure constant in QCD sum rule framework were already performed in 



a number of papers fT^ . The distinctive feature of the present work is a complete next-to-leading order 
analysis, containing correct treatment of Coulomb corrections. For discussion of importance of such 
corrections see, for example 

In QCD sum rules framework, the S^-meson structure constant is naturally obtained from two-point 
correlator of the following form 

n(g2) = z J d'xe''^-^{0\JB^{x)JBM^\0), (2) 

where Jb^ = bi'j^c. 

The left hand side of Eq.(2) can be computed in QCD, for much larger than Ag^^ (or , what 
alternatively, for — (m^ -|- much larger than A.q^^), with the use of short-distance Operator 

Product Expansion (OPE) for correlation function under consideration. 

Uiq%cD = npert{q') + Ca^iq'){-G') + ... , (3) 

71 

where Cc^iq"^) is a Wilson coefficient of gluon condensate operator and dots in right hand side of Eq.(|^) 
present contributions of operators with higher dimension [d > 4). 

The connection to physical spectrum of i?c-meson can be obtained by writing the following dispersion 
relation 

n(g2) = 1 / ^(flz^ + subtractions, (4) 
TT J s — q^ 



where 



P{s)had = TT— ^5(S - Ml) + p{s)pertO{s - Sthr), (5) 

[mb + rric) 

Here we have taken into account only lowest lying i?c-meson state, Sthr is a continuum threshold and 
p{s)pert is connected to Upertiq"^) via the following dispersion relation 

Upertiq'^) = — f ^^^^^'^^^ ^ + subtractions (6) 
TT J s — q^ 

In the numerical analysis we however explore a different anzaz, on which we will comment in the section 
with numerical results. 

There are several schemes of QCD sum rules, which can be used for the extraction of quantities, you 
are interesting in. The most popular among them are momentum and Borel schemes. In this paper we 
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will employ the first one and the studied object will be the momentum of two-point correlation function, 
given by the following expression 

^n = n!(^A^ n(g2)|,.=o (7) 

Thus far we have discussed QCD sum rule framework for determination of i?c-meson structure 
constant. As a primary goal of this paper is to perform a consistent analysis of the same quantity in 
NRQCD approximation, in the next section we define NRQCD sum rule framework and comment on 
connection of the latter with QCD sum rule analysis. 



3 NRQCD approximation 



In this section we set up a consistent NRQCD framework, in which the two-point correlation function 
n(g^) can be determined in a systematic manner at next-to-leading order. Our presentation in this and 
the next sections closely follow that of so we advice the reader to read that paper for more detail. 

NRQCD is an effective field theory of QCD designed to handle nonrelativistic heavy-quark-antiquark 
systems to in principle arbitrary precision. Considering all quarks of the first and second generation as 
massless the NRQCD Lagrangian reads [llH] 



•-NRQCD 
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bilinear terms and higher dimensional operators . 



The gluons and massless quarks are described by the conventional relativistic Lagrangian, where G^^, 
is the gluon field strength tensor, q the Dirac spinor of a massless quark and the gauge covariant 
derivative. For convenience, all color indices in Eq. (|^) and throughout this work are suppressed. 
The nonrelativistic c and b quarks are described by the Pauli spinors and respectively. Dt and 
D are the time and space components of the gauge covariant derivative D and = G*^* and = 
l^i-jkQjk electric and magnetic components of the gluon field strength tensor (in Coulomb gauge). 
The straightforward x^X bilinear terms are omitted and can be readily obtained. The short- distance 
coefficients ai, . . . ,05 are normalized to one at the Born level. The actual form of the higher order 
contributions to the short- distance coefficients dly . . . , IS irrelevant for this work, as we will later use 
the "direct matching" procedure at the level of the final result for the correlation function. 

Now let us discuss our correlation function in NRQCD approximation. Eq.@ can be rewritten as 

Uiq')=t{0\TJBMJBS-qV\0), (9) 
where JbX^) — Expressing Dirac fields for b and c - quarks in terms of Pauli spinors x ^"^^ 




(10) 
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we 



haveQ 



Y 8 V mfemc / / 

JsA-q)^ ^ -zf(^t^)(-g) - V^^^^^)'((DV^)t.Dx)(-g) + (13) 
Inserting these expansions into Eq. (^) we obtain 

^n(g^) = Ci(^hard, ^fact)-4l(-E', /^soft, /^fact) 
1 /TTif, — m, ^ ^ 



4 V mbrric 



C2(/ihard, /ifact)-4.2(i?, /igoft, /ifact) + . • • , (14) 



where 



A, = (0|(xV)(V^^x)|0), (15) 
A = ^(0|(xV)((D^)t-Dx)+h.c.|0), (16) 



The right-hand side of Eq. (p!4D just represents an apphcation of the factorization formahsm proposed 
in |]T5[. The second term in this expression is suppressed by v"^, i.e. of next-to-next-to-leading order and 
thus of no relevance to us in present analysis. The nonrelativistic current correlator Ai contains the 
resummation of the singular Coulomb terms. It incorporates all the long-distance, dynamics governed 
by soft scales like the relative three momentum ~ rriredV or the binding energy of the cb system ~ rriredV^- 
The constant Ci (it is normalized to one at the Born level), on the other hand, describes short-distance 
effects involving hard scales of the order of heavy quark mass. It is represented only by a simple 
power series in as and does not contain any resummations in a^. We would also like to note, that Ci is 
independent of q^. In Eq. (|l^) we have also indicated the dependence of the NRQCD correlators and the 
short-distance coefficients on the various renormalization scales: The factorization scale /ifact essentially 
represents the boundary between hard and soft momenta. The dependence on the factorization scale 
becomes explicit because of ultraviolet (UV) divergences contained in NRQCD. Because, as in any 
effective field theory, this boundary is not defined unambiguously, both the correlators and the short- 
distance coefficients in general depend on /Xfact- The soft scale /isoft and the hard scale /Xhard, on the 
other hand, are inherent to the correlators and the short- distance constants, respectively, governing their 
perturbative expansion. If we would have all orders in as and v at hand, the dependence of correlation 
function 11 (g^) on variations of each the three scales would vanish exactly. Unfortunately, we only 
perform the calculation up to next-to-leading order in as and v which leads to a residual dependence^ 
on the three scales /ifact, /isoft and /ihard- 



4 Calculation of NRQCD correlator 

To calculate the NRQCD correlator Ai we use methods originally developed for QED bound state 
calculations in the framework of NRQED [|l^, 0, |l^, ^ and apply them to Sc-meson bound state 

^Hcre and later in the paper, except stated otherwise, by nic and mf, we mean the pole heavy quark masses. 
^Here we would like to note, that the object studied P„ at NLO does not depend on fifact 
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description in the kinematic regime close to the threshold. At next-to-leading order quarks, composing 
Bc-Taeson experience only instantaneous interactions, given by the following potentials 



where 



T/(o) 

^ Utt 



ai 
rii 



[2(3o HM + ai 



— Ca Tnu 

3 3 



(17) 
(18) 
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Tni, 



(19) 



Here r = |r|, Cp 



3' 



3, T 



2' 



a. 



dsifJ'soit) and 7e is the Euler-Mascheroni constant. 



Thus, we can conclude that the problem of 5c-meson description close to threshold can be treated 
as a pure quantum two-body problem, so that we can use the well known analytic solutions of the non- 

i] and Rayleigh-Schrodinger time-independent 



relativistic Coulomb problem for positronium [20, 21 



perturbation theory (TIPT) to determine the corrections caused by all higher order interactions and 
effects. 

The calulational procedure for two-point NRQCD correlation function may be devided in the fol- 
lowing two steps 

Step 1: Solution of the Schrddinger equation. - The Green function of the next-to lead- 
ing Schrodinger equation is calculated incorporating the potentials displayed above. 
The correlator Ai is directly related to the zero-distance Green function of the 
Schrodinger equation. 

Step 2: Matching calculation. - The short- distance constant Ci is determined at 0{as) by 
matching QCD current Jb^ to corresponding NRQCD one. 

4.1 Solution of the Schrodinger equation 

The nonrelativistic correlator Ai is calculated by determining the Green function of the Schrodinger 
equation {E = \fq^ — {nih + mc)) 



+ 



\/i°)(f) + v;«(r) 



The relation between the correlator Ai and Green function reads 



A, 



6 



lim G{f, r', E) 



(20) 



:2i^ 



Eq. (^) can be quickly derived from the facts that G{f, r*, E) describes the propagation of b and c quark 
pair, which is produced and annihilated at relative distances |r| and 1^1, respectively, and that the same 
quark pair is produced and annihilated through the Jb^ current at zero distances. Therefore Ai must 
be proportional to lim|^j^/|^o G{f,f',E). The correct proportionality constant can then be determined 
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by matching the result of full QCD for perturbative spectral density , presented in Appendix A, to the 
imaginary part of the Green function of the Coulomb nonrelativistic Schrodinger equation. We would 



like to emphasize that the zero-distance Green function on the right hand side of Eqs. (|2T|) contains UV 
divergences which have to regularized. In the actual calculations we impose the explicit short-distance 
cutoff /ifact- As mentioned before, this is the reason why the correlators and the short-distance constants 
depend explicitly on the (factorization) scale /ifact- In this work we solve equation (^) perturbatively 
by starting from well known Green function Gf"* of the nonrelativistic Coulomb problem pO, pTl E2[ 



rrired 

and incorporate all the higher order terms using TIPT. 

The most general form of the Coulomb Green function reads (r = |r|, r' = |F|) 



(22) 



ATTT{l + ip)T{l-ip) 



1 oo 

dt / ds 



s{l-t) t(s-l) 



X 



1 



X 



dt ds 



ts 



sf—tf' 



exp I ip (^ \ \ {1 — t) + \ f\ {s — 1) + \ s f — t', 



where 



P = rriredV = JrrirediE + te 



P = 



2v 



r' < r, (23) 



(24) 



and r is the gamma function. The case r < r' is obtained by interchanging r and r'. Gf\f^ r*, E) rep- 
resents the analytical expression for the sum of ladder diagrams with Coulomb exchanges. The analytic 
form of the Coulomb Green function shown in Eq. ( 



) has been taken from Ref. |]20[. Fortunately we 
do not need the Coulomb Green function in its most general form but only its S-wave component with 
one of the relative quark distances set to zero.0 



Gf(0,r,E) = Gf\Q,r,E) 



27r 



dt e 



2iprt 



l+t 



— I 



2tt 



e'P' r(l - i p) U{1 -ip,2,-2ipr) 



(25) 



where U{a, b, z) is a confluent hypergeometric function |2^, It is an important fact that G^°''(0, r, E) 
diverges for the limit r because it contains power (oc 1/r) and logarithmic (oc Inr) divergences [p^ . 
As we have explained before these ultraviolet (UV) divergences are regularized by imposing the small 
distance cutoff pia.ct- The regularized form of limr^o Gc(0, r, E) reads 



G(°)'''^5(0,0,E) 



m 



red 



An 



iv — Cf Cb-i 



ln(— i 



/^fact 



2v 



(26) 



where the superscript "reg" indicates the cutoff regularization and '^{z) = dlnr{z)/dz is the digamma 
function. For the regularization we use the convention where all power divergences oc /ifact are freely 

^In the section, discussing the nonpcrturbative corrections, coming from gluon condensate operator, one more repre- 
sentation of Coulomb Green function wiU be introduced. 
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dropped and only logarithmic divergences oc ln(/ifact/?^red) are kept. Further, we define /ifact such 
that in the expression between the brackets all constants except the Euler-Mascheroni constant 7e are 
absorbed. The results for any other regularization scheme which suppressed power divergences (like the 
MS scheme) can be obtained by redefinition of the factorization scale. For convenience we suppress the 
superscript "reg" from now in this work. 

The Coulomb Green function contains cb bound state poles at the energies ^/s^ = nib + — 
Cpalmred/4:n'^ (n = 1, 2, . . . oo). These poles come from the digamma function in Eq. ( p6D and corre- 
spond to the nonrelativistic positronium state poles known from QED [^]. They are located entirely 
below the threshold point a/Sji^^ = m^+mc. This can be seen explicitly from the expression for imaginary 
part of Coulomb Green function G^\0, 0, E) 



Im 



n=l 

1 2 Cp O^gT^ 



where |\l/„(0)p = {mredCFds)^ /^T^n^ is the modulus squared of the LO nonrelativistic bound state wave 
functions for the radial quantum number n. The continuum contribution on the right-hand side of 
Eq. (^) is sometimes called "Sommerfeld factor" or "Fermi factor" in the literature. And the second 
term from the first line in Eq. ( ^Tf ) describes the resonance contributions. And finally, the corrections 
to the zero-distance Coulomb Green function calculated below lead to higher order contributions to the 
bound state energy levels poles and the continuum. 

Let us now come to the determination of the corrections to the zero-distance Coulomb Green function 
coming from the remaining term in the Schrodinger equation (pO]) . At next-to- leading order only the 
one-loop contributions to the Coulomb potential, V^^'^ have to considered. Using first order TIPT in 
configuration space representation the NLO corrections to Gf\Q^Q^E) reads 

G«(0,0,E) = - /rf¥G(°)(0,r,E)V«(r)Gf(r,0,E). (28) 



We will not calculate explicitly NLO correction to the Coulomb Green function here as the goal of 
this paper is to calculate the theoretical expressions for moments. The later can be most conveniently 
calculated by dispersion integration, using the following representation for the theoretical moments 

Pn=- -Im{Ci(^wd,/^so/i)a(0,0,E)}, (29) 

TT JO S" 

where E = ^/s — mi, — rric- 



4.2 Determination of the short distance coefficients 



The short-distance coefficient Ci and C2 can be determined by matching perturbative calculations of 
the matrix elements in full QCD and NRQCD. A convenient choice for matching is the matrix element 
between the vacuum and the state \cb) consisting of a c and a 6 on their perturbative mass shells with 
nonrelativistic four-momenta p and p' in the center of momentum frame: p + p' = 0. The matching 
condition is 



(0|67%5c|c6) 



{o\xli^c\cb) 



QCD 



+ C2 (0|(Dx6)^-DV^,|c6) 



NRQCD 



+ 



(30) 



NRQCD 
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To determine the short distance coefficients to order a^, we must calculate the matrix elements on 
both sides of ( |5UD to order a^. It is sufficient to calculate the order-a^ correction only for the coefficient 
Ci, since as we already said the contribution proportional to C2 is suppressed by v'^. The coefficient Ci 
can be isolated by taking the limit p — > 0, in which case the matrix element of {DxbV ' D'i/'c vanishes. 
Such calculations were done in 



28], where the following result was obtained for Ci. 



1 + 



: log 2 

nil) + nic nic 



(31) 



where mj-ed is the scale of the running coupling constant. The same result gives the direct matching of 
NRQCD correlator with QCD one, taking into account factor 2 for correction. 



5 Dispersion integration 

In this section we will discuss issues related to dispersion integration in expressions for NRQCD mo- 
ments. In general, the integration of spectral density over complete covariant form of integration 
measure is quite cumbersome However, in NRQCD approximation, as we will see soon, this task 
significantly simplifies. Let us make a change of variables E = ^/s — mi, — nic and consider a limit 
E <^ nih + rric. In this case the integration measure takes the form 

ds 1 2dE , , , , E 

exp{-(2n + l)ln(l + )} (32) 



{nib + nic)'^"' nib + nic nib + nic 

1 2dE , 2En , 2E , 
exp{ ■ } + 0{ ■ ) 



{mb + nic)'^"' nib + nic nib + nic nib + m 

The dispersion integration for NRQCD moments in this limit is 



P^!" = , ^ / -^eM-^^}KLoiE), (33) 

[nib + mcY"^ J E^„,i nib + nic nib + nic 

where Eund = "^''"''^^"^ is the negative binding energy of the lowest lying resonance. This integration 
is performed most efficiently by deforming the path of integration into negative complex energy plane, 
such that the part of the integration path parallel to imaginary axis is far away from bound state poles. 
That is 

= I ^ .,^ 7p / ^ exp{ }C,A,{E), (34) 

[mb + mc) 2t Jj-ioo nib + nic nib + nic 

where 7 <C Ebmd- Note, that in the above equation the real part of correlator Ai is also present, which 
is needed for integration over the new path. After performing the second change of variables E —E 

,^ TT 1 r'+'°° 2dE , 2En 



/ ^^e^v{^^}C,A,{-E) (35) 
Ji-ioQ nih + nir nib + nir 



{nib + nicY"^ 2TTi J-y-ioo nib + nic nib + 

we see, that the problem of evaluation of NRQCD moments can be related to the inverse Laplace 
transform of integrand expression, for with there are a lot of tables in literature. The procedure of 
taking inverse Laplace transform can be further simplified by noting that integration path is far away 
from bound state energies and hence the integrand can be safely expanded in as- 
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Following the steps, described above, and using relations for inverse Laplace transform from Ap- 
pendix B the NRQCD moments in leading order of NRQCD expansion have the form 



where 



and 



2ak;;.)^^;3/. ^"(7), (36) 

$"(7) = l + 2v/(vr)7 + ^7' + 4v^f:Q exp{ Q }[l + erff^)] (37) 



{nib + rricy/^ 



The calculation of NLO order correction to moments, coming from correction to potential (|TB[) is a bit 
more involved. With the help of presentation (^) for Coulomb Green function we have 

^j^thiNLo 6 1 n+^oo 2dE r ^ 



{rrib + mcY'^ 2i J^-ioo mi, + rric mb + rric 

2 



Jo \ TV J Jo Jo \ tU 




(39) 



where A = "Ved"" and k = 2mredE. The integration over r can be easily performed explicitly, while 
the situation with integrations over t and u is far more complicated. However, as we noted above the 
integrand expression in the case under consideration can be expanded in series over A| . Then the result 
of integration over all terms in A expansion, except A° can easily expressed in terms of functions Wp 
and Wp, introduced in [^0. As for the term with A°, which contain a manifestly divergent integral, 
it can be easily calculated in terms of K{t) = (0| exp(— i7r)|0), using the free evolution function 
Ko{t) = {0\exp{-HoT)\0). Here H = Hq - Hq = and r is an Euclidean time. The 

evolution function Kir^ can be related to the n'th NRQCD moment by the following relation 



= 7 ' K{^^) (40) 

Using this relation and explicit expression for free quark propagation function: 

(x|exp(-ifor)|y) = (^)'^'exp (-!^(x-y)^) , (41) 

we can now easily evaluate the contribution of A° term for NRQCD moments. The combined result for 
NLO order correction to NRQCD moments reads 

tpthlNLO ^ 3(mr.ed)^^^(mb + "^c)^^^ ^W N (Ar)\ 

^ 2v/?(m6 + me)2"n3/2 ^^^^ 



^It's the same as the expansion in as 

®The expressions for them can be found in Appendix B 
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where 

Here Hi = jj-aoft exp (^) • The full perturbative result for NRQCD moments with account of hard gluon 
corrections reads 



j-pi/ijpert — _|_ '^sy^rcd/ 



mb-rric , nib 

log 2 



)^[pt.]LO^[p*.]A.LO| (44) 



6 Nonperturbative corrections 

In this subsection we would like to consider corrections, given by gluon condesate operator. The 
calculation is very similar to the one done previously by Voloshin and Leutwyler for the case of 
equal quark masses. For the determination of corrections to Coulomb Green function, coming from 
gluon condensate operator, we will exploit the fact, that the size of vacuum fluctuations of gluon field is 
much larger then the size of i?c-nieson state. So, we can perform a multipole expansion for interaction 
of i?c-nieson with gluon condensate, whose first term is 

^int = -IgCrE^, (45) 

where f = Xc — xi, = Anas, ^'^ = — ^2- emploing colour and Lorenz invariance of the vacuum 
state one can relate the average over the vacuum state of two chromoelectric fields to the manifestly 
Lorenz-invariant value 

1 

96' 

Thus, the Coulomb Green function with account for gluon condesate corrections has the following form 



\EtEl\0) = -_5,,5''^(0|G^,G^J0) (46) 



Gix,y,E) = G(o)-^(0|7ra,G;:,G;jO)x (47) 



dVdV(rV)G(o)(f , r-, E)G^s)ir, r"^, E)G^o)ir, y, E), 

where G(o) and (7(8) are Coulomb Green functions in singlet and octet colour states correspondingly. 
They are defined as solutions to the following equation 



where 

The expressions to NRQCD moments, coming form gluon condensate, can be calculated with the 
use of wave decomposition of NRQCD Green function 

00 

Gf^'\x,y-E)=Y.i'^l + l)G^^''^\x,y-E)mxy)l^y), (50) 
z=o 
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where 



TT ' ' ^ ^t^o(s + / + l-m,eda(0'8)/A;)(s + / + l)!' 

(51) 

where x = \x\, y = \y\. Pi and are Legendre and Laguerre polynomials. The result for gluon 
condensate correction to Coulomb Green function is[]: 

[Pn\ - i2(m, + m,)2n+3 x(7)(0|a.G^.G^J0) (52) 

The function x(7) is given by the following equation 

. 4V?{|:(p+l)(p + 2)(p + 3){^fe(2)-^,(2,(l + ^(25-^))l 

^ " ■)(l + ^(5-;^))l+ (53) 



9p + 17'^ >+l' >+l'' 6 ' 9n + 17' 



p + 2'" > + 2' 18" > + 2'' 9p+19'" > + 3' > + 3' 

^(5 + ^^))] + ^—[U^) - M^){1 - ^(25 + 
6 ^ 9p+ 19^^^ 9p + 20^^ > + 4^ ^^> + 4^^ 12 ^ 

6 4 °° — 1 7 

9j9 + 20^^^^ + ^ E (81p2_ i)(8ip2_4)<^i(-^)}' 



;.,(x) = a;-3[e^''(i + erf(x))-l]-(4=)a;"'-a;-\ (54) 



where 



02(x) = [e^ (1 + erf(x)) - Ijx"^ - ^x). (55) 

The above expression in the limit of equal quark masses coincides can be checked with the one derived 
previously by M.B.Voloshin However, this expression is quite complicated and in the range 7 < 1.5 
it is far more convenient to use an approximated formulum x{l) = e~°'^°'''$°(7). 



7 Numerical results 

In this section we will discuss our numerical estimates for various c-quark mass definitions, the i?c-meson 
mass and coupling constant. Let us first discuss our anzaz for experimental spectral density, the need 
for which is dictated by our present lack of experimental measurements on higher i?c-meson states. The 
experience in potential models and quasi-local sum rules |^T| allow us to write 



M„ = Ml + 2rinn, Mi = 6.3 GeV, T = 0.415 GeV, 



2 - ^ , /i(0-) = 0.4GeV, (56) 



^Here we assume that one should perform the dispertion integration in order to obtain the gluon condensate correction 
to NRQCD moments 
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where M„ is the mass of nS i?c-nieson state, and is its couphng. So that, the experimental spectral 
density has the form 

3 r2nj4 

P{s)hard = 7^Y1 r^^ir^^(" - ^n) + Pis)pert0{s - S*^.) (57) 
n=l '^cj 

Now with the knowledge of experimental spectral density and thus of experimental moments we can 
compare them with theoretical ones. To begin with let's see what values of the pole c-quark mass will 
give us the required values of experimental moments. For the value of 6-quark pole mass M^"'^, needed 
for calculation of theoretical moments^, we have taken the result of NNLO analysis for T family ||Tl|, so 
that we have Mf"'^ = 4.8 ±0.06 GeV. Varying the soft scale in psoft the range 1.1 — 1.2 GeV and fixing 
the hard scale Phard to 2. GeV, from the stability of the ratios of experimental to theoretical moments we 
have that M^"'*^ = 1.7 — 2.1 GeV. The range of pole c-quark mass obtained is quite large and to reduce 
it we need some extra condition to satisfy. Again, from the spectroscopy of i?c-meson family it is known 
that the characteristic value of strong coupling constant for heavy quark dynamics inside Sc-meson is 
in the range Os = 0.43 — 0.48. This extra requirement heavily constrains the value of c-quark pole mass, 
so that M^"'^ = 2.03 ± 0.06 GeV. At Fig. |I| we have shown the ratios of experimental to theoretical 
moments as functions of momentum number at fixed value of Hsoft = 1-1 GeV and different values of 
c-quark pole mass. 




[ ] 

5 10 15 20 25 n 

Figure 1: Ratios of experimental and theoretically calculated moments -jjur as functions of momentum 
number (The solid line is for M^°'^ = 2.03 GeV, the line with short dashes is for M^°'^ = 1.9 GeV and 
the one with long dashes is for M^"'*^ = 2.1 GeV. All ratios are plotted for fihard = 2.0 GeV and 
^isoft = 1.1 GeV ). 



®The precise values of other parameters, that is gluon condensate and continuum threshold are not so important, as 
the correlation function depends on them weakly. 

^Note, that at NLO, due to vanishing of anomalous dimension of the current under consideration, the dependence of 
two-point correlator from fihard enters only through as{fJ-hard)- So, the particular value of ^uard is not very important. 
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The performed analysis shows that to satisfy chosen criteria we need very low values of fisoft ~ 1.1 GeV 
or lower. On the other hand, as was shown in to match perturbative heavy quark potential with 
full QCD potential at characteristic distances of quarks bound inside heavy-heavy mesons, we should 
have the soft scale fisoft in the region 1.3 — 2.0 GeV[3 To achieve this goal we modify the original TIPT 
for theoretical moments by shifting the leading order Coulomb potential by a part of constant term in 
NLO potential, that is 



(0) 



modified 



47r 



(5^ 



The resulting theoretical moments will now depend not only from hard and soft scales but also from 
the value of the shift parameter t. To get rid of the latter dependence, while keeping in mind our goal, 
we perform optimization in parameter t, require the NLO corrections to theoretical moments at given 
momentum number n do not exceed 1%Q. As a result we get M^"'^ = 1.96 ± 0.05 GeV. At Fig. |T] 
we show the ratios of experimental to theoretical moments as functions of momentum number at fixed 
value of ^soft = 1-3 GeV and different values of c-quark pole mass. We see, that within the error bars 
our results for the c-quark pole mass from these two estimates agree with each other as well as with the 



results of other estimates 
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Figure 2: Ratios of experimental and theoretically calculated moments -thtt as functions of momentum 
number in modified TITP (The solid line is for M^"'^ = 1.96 GeV, the line with short dashes is for 
][{poie _ j^ gg GeV and the one with long dashes is for M^"'"^ = 2.05 GeV. All ratios are plotted for 
l^'hard = 2.0 GeV and fi^oft = 1-3 GeV ). 



^"Such choice of soft scale is also desirable in order to have reliable perturbative predictions for theoretical moments. 
^^Note, that here we will need already ai,(^so/t) = 0.43 — 0.48 
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Having extracted the value of c-quark pole mass from Bc-snm rules, an analogous analysis can be 
performed for other definitions of heavy quark masses. Below we give the numerical values of these 
masses obtained in the same way and comment on their relation with each other. First, let us consider 
the running c-quark mass related to the pole one by the following relation ^ 

1 + 1,333 f^U 113^44 -L041,,/"-''"">^' 



171(171} \ TT / \ VT 



[190.1 - 26.7ni + O.QbSnJ] ( ^iM^ (59) 



TT 



Our estimates for the m{m) mass are 

me(mc) = 1.40 ± 0.07 GeV (60) 

However, not pole not running quark masses are qood ones when one works in the region of energies close 
to threshold. As we have seen in the analysis performed above the two-point correlation function has 
a strong dependence on renormalization scales as well as correlation between the pole mass and strong 
coupling constant (the situation is similar for running quark masses), so to reduce this dependence 
somewhat it's would be more appropriate to explore quark masses, whose use can eliminate as much 
as possible all such dependencies and correlations among parameters. For this reason here we give 
estimates of IS* and potential subtracted c-quark masses. The former can be related to the pole mass 
jYfpoie -^[^i^ ^i^Q help of following formulae 



where 



= (62) 
= (f)[/5o(^ + l) + y], (63) 



,2 



'^2 , 2i 



16 8 ^ 48^ ^ V ; 

L = ln(- — ^ r), (65) 



and 



/4343 . vr^ 22 \ ^2 /1798 56 

aa = V^Ti \ Cs \C\-\ \ Ca CaTui 

V 162 4 3 ^7 ^ V 81 3 

^±-lQC,)CpTni+(^Tn^' (66) 

With the use of these formulae we have the following estimate for IS* c-quark mass at optimized value 
of /i.o/t = 1.3 GeV|3 

M^^ = 1.49 ± 0.07 GeV (67) 

^^This value is taken from the previous analysis for c-quark pole mass performed in modified TIPT, as the shift of 
leading order Coulomb potential implicitly includes part of NNLO corrections. However, as we do not have up to date 
the complete NNLO order analysis this result should be taken with care. 
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Having done this estimate it is instructive to compare it with the estimates of NLO analysis done in 
IS scheme. To obtain the expressions for theoretical moments in this scheme we make a following 
substitution 



1 



1 



exp-2"^^°("^){l-2n 



2n 



(6^ 



In all other places the heavy quark pole masses should be changed to IS* masses. Varying the scales in 
the same ranges as we did for pole c-quark mass we obtain, that the sum rules are stable themselves, but 
the ratio of experimental to theoretical moments is far from unity for reasonable values of IS" c-quark 
massQ. To solve this problem we again consider the modified TIPT for theoretical moments. It is 
important, when writing expressions for theoretical moments, to make for consistency a similar shift in 
NLO order relation between pole and IS* masses. Performing calculations in modified 15* scheme we 
obtain the following estimate for 15* c-quark mass: M]^ = 1.52 ± 0.05 GeV, which, within the errors, 
is in agreement with previous estimate for 15* c-quark mass. Fig. ^ shows us the ratios of experimental 
to theoretical moments as functions of momentum number at fixed values of /iso/t = 1-3 GeV and 15 
c-quark mass. 



1 . 2 r 



1.15- 



1.1- 



1.05- 




35 n 



Figure 3: Ratio of experimental and theoretically calculated moments -jfth- as functions of momentum 
number in modified TITP {M^^ = 1.52 GeV, fihard = 2.0 GeV and fisoft = 1-3 GeV ). 



Recall that this mass by definition differs from a half of J/\E' mass on a small value given by nonperturbative 
corrections. 
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To finish the discussion of 15* c-quark mass we note, that even this mass was defined as M^"'^ minus 
one half of perturbative Coulomb energy in J/\E'-meson, it can be equally well applied for i?c-nieson 
as well as other systems, containing c-quark. The arguments here are the analog of Upsilon expansion 
for charmed quarks and the fact, that the Coulomb energy for static quarks does not depend on their 
flavors. 

The estimate for the value of potential subtracted c-quark mass was obtained from the relation of 
the latter to the running quark mass 



mpsifJ^f) = m(m){l + 



■[1- 



+ 



m[m) 



TT 



[13.44- 1.04nz 



^(Gi + 47r/3o[ln 



3m{m 



e\2 



2])]} 



Thus, at /i/ = 1.5 GeV we have 



mp5(1.5GeV) = 1.42 ± 0.07GeV 



(69) 



(70) 



Having completed the estimates of different c-quark masses, we can now perform estimates of the mass 
and coupling constant for the ground i?c-nieson state, fixing the heavy quark masses at their central 
values. 



6. 304 



6.303 - 



6 . 302 - 



6 . 301 - 




Figure 4: The values of Be - meson mass extracted from the two-point NRQCD sum rules in moment 
scheme. 



The results for mass and coupling of Sc-meson in momentum scheme of NRQCD sum rules can be 
easily seen from Fig. § and |^. 
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Figure 5: The values of - meson coupling constant extracted from the two-point NRQCD sum rules 
in moment scheme. 



Recently, the estimate of Be mass was obtained in the perturbative potential approach by fitting the 



masses of J/^ and T in order to get a good covergency in as and extract the heavy quark masses [ p^ . 
The perturbative mass = 6.323 ± 0.007 GeV is very close to the our estimate in the framework 
of QCD sum rules. This fact indicates a small nonperturbative correction, which was limited by —60 



MeV in p6| 



8 Conclusion 

In this paper we have presented complete NLO analysis for i?c-meson two-point NRQCD sum rules. 
Analitycal results for perturbative spectral density and gluon condensate contribution with account for 
summed Coulomb corrections are derived and analyzed. A detail numerical analysis as well as discussion 
on the determination of various c-quark masses together with couplings and masses of lowest lying Bc- 
meson resonances from the mentioned sum rules are provided. The analysis shows that to reduce 
the uncertainties of calculations it is mandatory to have complete NNLO expressions for theoretical 
moments, which we plan to accomplish in nearest future. 
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Appendix A 

In this Appendix we present the two-loop result for perturbative spectral density of Be meson P3| : 




log ai log A 



(71) 



i=b,c 




where 




(72) 



and 




19 + 2^2 + 3f ^ 
32^^ 



mi{mi — rrij) 



q^v{l + v) 
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1+V + 
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) 



(73) 
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Vl + at 




(74) 
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Appendix B 



In this Appendix we have collected some formulae, needed for calculation of NRQCD moments in next 
to leading order. 



1 /■7+ioo 1 
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r(i) 



dtl c^«l^Mi±^ln- 

^(2 + m) " 



tu J 

'{i + t){i + uy 

tu 



(p + i)Cp+i 
r(|) ^ 



m] 



where is the Rieman zeta function for argument p. 
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